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Abstract 

The  paper  is  devoted  to  the  study  of  the  long-4erm  stability  of  the  propagation  delay  of  a  fiber-- 
optic  timing  link  under  real  environmental  conditions.  The  temperature  dependence  of  the  light 
propagation  delay  on  surrounding  temperature  changes  has  been  found  as  dominant  An 
appropriate  method  has  been  developed  to  correct  mathematically  the  thermal  delay  variations 
in  the  fiber.  After  the  correction  for  this  temperature  effect^  only  subnanosecond  residual 
variations  have  remained. 

Introduction 

To  ensure  subnanosecond  time  comparisons  of  the  remote  clocks  over  a  long  period  of  time  (one  year  and 
more),  an  adequate  transmission  system  must  be  used  in  addition  to  stable  measurement  and  distribution 
systems.  A  time  transfer  link  based  on  the  optic  fiber  seems  to  fulfill  the  requirements  for  precise  timing 
over  a  distance  of  several  tens  of  kilometers.  The  electronic  part  of  the  optic  transfer  link  can  be  basically 
well  designed  and  maintained  in  closely  controlled  laboratory  conditions.  Also,  the  parameters  of  the 
electronics  can  be  regularly  checked  to  obtain  good  long-term  time  transfer  stability.  In  contrast,  the 
fiber,  which  spans  the  outer  distances  between  the  laboratories,  is  exposed  to  uncontrollable  natural 
environment  even  if  buried  under  ground,  taking  into  account  the  insensitivity  of  the  fiber  to  the 
electromagnetic  perturbation,  all  other  environmental  factors  like  temper^ure,  humidity,  atmospheric 
pressure,  etc.  may  have  influence  upon  it.  All  these  factors  may  cause  the  instability  of  the  propagation 
delay  of  the  optical  signal  in  the  fiber.  The  spurious  fluctuations  of  the  link  propagation  delay 
superimposed  on  true  measured  phase-time  values  could  be  mistakenly  interpreted  as  the  instability  of 
the  compared  clocks  unless  specif  care  (correction,  compensation)  is  taken  with  respect  to  the  measured 
data. 

For  this  reason  the  investigation  of  time  transfer  parameters  of  fiber-optic  links  was  one  of  the  first  tasks 
which  arose  just  after  the  realization  of  the  Czech  National  Time  and  Frequency  Group  Standard 
(CNTFGS)  [1]  in  1993.  The  CNTFGS  is  based  on  the  HP  5071A  cesium  clocks  operating  in  two  remote 
laboratories  in  Prague.  Fiber-optic  links  are  used  to  provide  the  internal  traceability  of  individual  clocks. 
These  transfer  links  (  each  of  than  consisting  of  an  idaitical  1300  nm  optical  transmitter,  approximately 
10-km4ong  single-mode  fiber  buried  underground,  and  an  optical  receiver  followed  by  a  5  MHz  PLL 
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filter)  are  used  to  transfer  5  MHz  ou^ut  signals  of  the  clocks  in  both  directions.  The  two-way 
configuration  of  the  CNTFGS  time  transfer  network  permits  use  of  the  standard  data  obtained  from  the 
routine  time  scale  comparisons  for  evaluation  of  the  link  transfer  stabihty.  This  paper  summarizes  the 
results  of  nearly  four-year  study  of  the  instability  of  the  propagation  delay  of  the  fiber  timing  links.  It 
has  been  found  that  the  temperature  of  the  soil  surrounding  the  optic  cable,  which  causes  the  tharmai 
dilatation  of  the  fiber,  is  dominant  among  all  other  environmental  factors  which  also  might  affect  the 
long-term  propagation  delay  fluctuations  of  the  whole  link. 

The  experiment  configuration  and  time  relations 

The  basic  configuration  of  the  clocks  and  timing  links  is  depicted  in  Figure  1.  One  cesium  clock  is  located 
in  laboratory  A  and  two  clocks  in  laboratory  B.  Identical  automated  measuring  systems  are  used  at  both 
laboratories  to  provide  simultaneous  phase-time  measurements.  Two  output  signals  transmitted  from  B 
at  the  time  instants  Ti  and  T2  using  the  links  1  and  2  are  received  at  A  at  the  time  instants  Ti*  and  T2*, 
respectively.  In  the  same  way,  the  output  signal  transmitted  from  A  at  the  time  instant  To  using  the  link  0 
is  received  at  B  at  To*.  Each  transfer  link  is  characterized  by  its  total  propagation  delay  which  can  be 
modeled  as  a  sum  of  a  constant  and  a  variable  parts.  Only  these  variable  components  t  i ,  (i  =  0,1,2)  will 
be  of  further  interest .  To  make  the  remote  clock  comparisons,  the  differences  To  -Ti  and  To  -T2  should  be 
known.  Apparently,  they  cannot  be  measured  directly,  but  they  can  be  evaluated  from  the  differences 
measured  at  laboratory  A: 


ATai^T,-Ti*  (la) 

A  Ta2  =  T,  -  T2*  (lb) 

ATa12=Ti*-T2*  (Ic) 

and  at  laboratory  B: 

ATbi  =  To*-Tj  (2a) 

ATb2  =  To*-T2  (2b) 

^Tbi2  —  Tj-  -  T2  ■  (2c) 

Since  T*  =  Ti  +  Ti,  the  desired  phase-time  differences 

To-Ti  =  ATai  +  r?  =  ATbi~  ro  (3a) 

T0-T2  =  ATa2  +  72  =  ATsi-ro  (3a) 

and  the  following  relation  can  be  derived  from  (1)  and  (2). 

ATbi  -  ATai  =  zlj  =  To  +  r?  (4a) 

A3'b2  -  ATa2  ~  A2=  To  +  t2  (4b) 

ATbi2  -  ATai2  ~  S  =  T2  -  Tz  (4c) 

ATai  +  ATbi  =  2(To -Tj)  +  (tq- Ti)  (5a) 

ATa2  +  a  Tb2  2(To  -TV  +  (to  -  T2).  (5b) 


Ai  and  A2  denote  the  summary  delay  of  the  loops  1  and  2  made  up  by  the  link  0  together  with  the  links  1 
and  2  respectively,  and  5  denotes  the  differential  delay  of  the  links  1  and  2. 
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Analysis  of  obtained  data 


DELAY  OF  THE  LOOPS  1  AND  2  AND  ITS  CORRECTION 


The  plots  of  daily  values  of  Ai  and  A2  are  displayed  inFigure  2.  One  can  see  that  the  curves  exhibit  very 
■similar  seasonal  variations  with,  the  peak-to-peak  value  being  approximately  15  ns.  The  maximum 
frequency  deviation  is  about  +  1.10"'^.  The  third  curve  displayed  inFigure  2  represents  mean  daily 
temperature  Tout  taken  in  Prague  .  The  close  correlation  of  the  long-term  propagation  delay  with  the 
outdoor  temperature  is  obvious.  The  detailed  study  shows  that  the  correlation  coefficient  reaches  its 
maximum  value  (  p  =  0,85)  if  the  Ai  or  A2  functions  are  delayed  by  about  20  days  with  respect  to  Tout  ■ 
Such  a  delay  of  the  response  to  the  outdoor  temperature  might  be  explained  only  by  the  slow 
propagation  of  heat  in  the  soil.  Therefore,it  can  be  assumed  that  the  underground  temperature  which  has 
an  influence  on  the  fiber  is  the  source  of  this  type  of  fluctuations.  To  confirm  this  assumption,  the 
underground  temperature  along  the  cable  should  be  known.  However,  the  direct  measurement  of  this 
temperature  is  alrnost  impossible.  The  fact  is  that  the  9-km-4ong  section  of  the  optic  cable  traces  the  urban 
area  and  is  buried  in  different  depths  varying  fi'om  0.5  m  to  3  m  (part  of  the  cable  is  led  in  a  cableduct ; 
another  part  is  laid  in  the  soil).  So  the  soil  temperature  was  modeled  under  the  following  simplifying 
presumptions: 

-  the  soil  was  presumed  to  be  homogeneous, 

-  the  mean  daily  temperature  of  the  soil  surface  was  supposed  to  be  equal  to  the  mean  temperature  of  the 
air, 

-  there  are  no  other  sources  of  heat  influencing  the  optic  cable. 

Under  these  conditions  the  soil  temperature  indifferent  depths  0.25  m  to  3  m  was  evaluated  as  described 
in  [3]  .  The  correlation  of  the  Ai  and  A2  with  the  relative  daily  values  of  soil  temperature  in  different 
depths  was  analyzed, yielding  the  results  displayed  inFigure  3.  It  is  seen  that  the  maximum  correlation 
coefficients  (  Pai.t®  =  0.93  for  Ai  and  pA2,T(h)  =  0.9  for  A2)  are  obtained  for  the  depth  h  =  1.3  m. 

An  attempt  has  been  made  to  find  a  correlation  of  the  A 1  and  A2  with  other  environmental  factors  (air 
pressure,  indoor  temperature  of  the  laboratories  A  and  B).  None  of  the  correlation  coefficients  have 
exceeded  0.05.  This  indicates  that  the  temperature  of  the  soil  surrounding  the  optic  cable  has  the 
dominating  influence  on  the  long-term  stability  of  its  propagation  delay  . 

In  order  to  make  the  optimal  correction  of  the  temperature  changes  of  the  link  delay,  the  estimation 
method  has  been  used.  The  temperature  coefficient  K;  [ns/  km.°C]  minimizing  the  value  A(h)  for  different 
values  of  the  depth  were  chosen  for  the  loop  1  and  2: 

av 

^i{n)-Ki.d.n{n)\  (6) 

where  A(i)  -samples  of  daily  values  of  the  loop  delay 

Th(n)  -  samples  of  daily  values  of  the  soil  temperature  for  the  depth  h 
Ni  =  MJD^ 

N2  =  MJDnmt . 


i  =  1, 2  denotes  the  number  of  loop 
d  -  length  in  kna  of  the  buried  part  of  the  cable  ( d=9  +  9=  18) 

The  minimum  A(h)  =  1  ns  was  found  for  h  =  1.3  m .  That  is  in  a  close  agreement  with  the  correlation  result 
mentioned  above.  The  values  Ki  =  0.046  andK2=  0.041  obtained  from  (6)  for  this  ^effectiv^  depth  do 
not  differ  very  much.  This  result  is  not  surprising  if  we  take  into  account  that  the  individual  fibers  placed 
within  one  cable  should  exhibit  very  similar  propagation  delay  wiations  under  the  same  environmental 
conditions.  That  makes  it  possible  to  simphty  the  above  relations  by  assuming 

To  =  Ti  =  Tj  =  T.  (7) 

Under  this  assumption 

A=Aj  =  A2  =  2  T  (8) 

holds  and  so  the  temperature-corrected  loop  delays  Aicar  (  i  =  1, 2 )  can  be  evaluated  as 

Aico.  =  Ai-K,d.n  (9) 

The  plots  of  corrected  time  fimctions  Aicm  and  Aaeon  are  shown  in  Figure  4  .  They  have  no  long-term 
variationSt  but  they  show  residual  noise  fluctuations  as  the  curves  Ai  an  A2  in  Figure  3.  Note  that  the 
records  inFigure  3  and  4  consist  of  two  different  parts.  The  left-hand  part  (MJD  49379  to  MJD  50299)  has 
been  obtained  from  the  data  measured  by  less  accurate  time  interval  counters,  which  provided  uncertainty 
of  CT  =  0.52  ns.  Since  tiiat  uncertainty  was  too  large,  high  resolution  counters  with  a  =  30  ps  were  used 
later  in  the  course  of  ejqjeriment  (MJD  50300  to  MJD  50687). 

Among  the  measures  used  in  the  time  domain,  the  time  difference 


WEV(T)  =  ^^((jfxf}  (10) 

seons  to  be  optimal  to  describe  the  time  stability  of  the  uncorrected  time  transfer  and  to  estimate  the 
effect  of  temperature  correction  .TDEV  (1  day  ^  x  <  240  days)  was  evaluated  first  for  A2  within  the 
time  interval  MJD  49379  to  50299.  The  corresponding  plot  daioted  as  A2  LRC  is  presented  inFigure  5. 
It  is  seen  that  the  TDEV(x)  is  basically  constant  for  1  day  <  x  <  30  days.  The  value  of  TDEV 
(x=lday)  =  0.94  ns  .  The  TDEV(x)  is  increasing  for  x  >30  days  and  it  reaches  its  maximum  for  x  =  180 
days.  This  is  the  result  of  the  seasonal  fluctuations.  Then  the  TDEV  (1<  t  <  90  days)  was  evaluated  for  A2 
within  the  time  interval  MJD  50300  to  MJD  50687.  In  this  case  the  curve  A2  HRC  representing  TDEV(x) 
characterizes  the  instability  of  the  loop  much  better  because  of  substantially  higher  resolution  of  the 
counters.  The  plot  A2HRC  starting  from  TDEV2  (1  day)  =  0.4  ns  is  monotonously  increasing  and  reaching 
the  A2LRCforx=  90  days. 

The  corresponding  TDEV(x)  plots  (denoted  by  A2cotr  LRC  and  A2c(ar  HRC)  calculated  for  A2coir  in  the 
above  mentioned  intervals  of  x  and  MJD  are  displayed  inFigure  6.  We  can  see  that  the  TDEV(x)  values 
characterizing  the  temperature-corrected  loop  propagation  delay  well  correspond  to  the  ones  of 
uncorrected  loop  propagation  delay  within  the  interval  1  day  ^  x  <  30  days  .  The  correction  effect  is 
evident  for  x  >  30  days. 
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CORRECTION  OF  MEASURED  PHASE  TIME  DATA 


Knowing  the  law  of  temperature  dependence  of  the  propagation  delay,  we  can  conect  for  the  measured 
phase-time  data  in  order  to  remove  the  long-term  fluctuations  due  to  the  transmission  link.  For  one-way 
time  transfer  the  following  formula  (i=l,2)  can  be  used: 

TrTi^ATAi-K.d.Ty,,  (11) 

This  correction  needs  the  current  soil  temperature  Th  to  be  evaluated  from  the  mean  daily  outdoor 
temperature,  as  described  in  [3] .  Let  us  mention  that  the  outdoor  temperature  changes  will  exert  an 
influence  upon  the  fiber  of  the  described  system  after  a  delay  of  20  days. 


Another  way  to  correct  the  data  is  to  take  advantage  of  the  two-way  transmission.  It  is  clear  (if  (7)  and 
(8)  are  valid)  that  the  desired  time  differences  cleared  of  transmission  link  fluctuations  can  be  obtained 
simply  by  using  (5a,  5b)  e.g. 


To-Ti  = 


ATai  4-  AThi 


(12) 


even  without  the  knowledge  of  the  actual  values  of  link  propagation  delay. 

Conclusion 

The  obtained  results  demonstrate  the  possibility  of  stabilizing  mathematically  the  propagation  delay  of  the 
studied  time  transfer  optic  link  within  1  ns  or  less  (in  the  sense  of  TDEV(t))  fortune  intervals -r  >1 
day  by  correcting  for  the  influence  of  the  soil  temperature  variations  upon  the  optic  fiber  buried  under¬ 
ground.  This  correction  reduces  the  long-term  instability  for  time  intervals  longer  than  30  days  and 
permits  us  to  reach  1  ns  timing  uncertainty  within  the  group  standard.  In  the  course  of  almost  four-year 
period  of  experiment,  no  significant  systematic  phase  drifts  (aging, etc.)  have  been  observed. 

Even  after  the  temperature  correction  residual  phase  fluctuations  for  shorter  time  intervals,  t  £  30  days, 
remain.  Their  source  is  not  yet  fully  known;  however,  from  the  data  ofthe  two  sets  of  counters  with 
different  resolution  we  can  conclude  that  at  least  a  part  of  the  fluctuations  is  due  to  the  limited  counter 
resolution. 
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Fig.  1  :  Basic  diagram  of  the  time  transfer  system.  One  dodc  signal  is  transmitted  from  the  laboratory  A 
to  the  laboratory  B  and  two  clock  dgnals  are  transmitted  in  a  opposite  direction  using  idrartical  fiber-optic 
transfer  links. 


Fig.  2  This  figure  shows  the  plots  of  relative  summary  propagation  delay  A]  of  the  loop  1  created  by 
Imks  0  and  1,  and  A2  of  the  loop  2  created  by  links  0  and  2.  The  third  plot  displays  the  mean  daily 
outdoor  temperature  Tout  in  Prague. 
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Fig.3  Correlation  of  loop  delays  with  imdergroimd  temperature  in  different  depths.  The  correlation 
coefficient  p  reaches  its  maximum  value  for  the  depth  h  =  1.3  m. 
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Fig.4  displays  the  plots  of  the  relative  propagation  delays  of  the  loops  1  and  2  after  correction  for  the 
soil  temperature  in  the  depth  h  -  1 .3  m 
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Fig.  6  displays  the  plots  of  TDEV  (t)  of  corrected  propagation  delay  of  the  loop  2  calculated  from 
the  data  measured  by  low  resolution  counters  (A2  con,  LRC)  and  by  high  resolution  counters 
(A2C017  HRC). 


Questions  and  Answers 


ALBERT  KIRK  (JPL):  You  show  a  seasonal  temperature  plot  and  you  also  mention  that  the  depth  of  the 
cable  was  between  .5  (point  5)  and  three  meters.  At  what  depth  did  you  measure  the  temperature? 

OTOKAR  BUZEK:  The  temperature  measurement  was  not  down  under  the  ground.  We  used,  as  input 
data,  the  mean  temperature  of  Ae  air  in  Prague.  It  was  supposed  that  this  is  the  ground  surface  temperature. 
Then  using  the  standard  equations,  the  temperature  imderground  was  calculated.  It  was  only  calculated, 
because  the  soil  is  operating  as  a  low-pass  filter.  Only  the  long  wave  of  the  heat  is  transported  down  and 
time  delayed.  So  it  was  not  measured,  it  was  calculated. 

ALBERT  KIRK:  Was  most  of  the  cable  at  the  three-meter  depth  or  between  .5  (point  5)  and  three  meters? 
Can  you  express  that  in  a  certain  percentage,  the  relationship  of  depth  for  the  10  kilometers? 

MALCOLM  CALHOUN  (JPL):  May  I  answer  that? 

ALBERT  KIRK:  Yes,  please. 

MALCOLM  CALHOUN;  Eight  hundred  meters  was  at  a  half-meter  depth.  The  rest  of  the  link  was  at  two 
to  three  meters.  Unfortunately,  the  800-meter  at  a  half-meter  depth  would  account  for  the  biggest  part  of 
the  cycling  that  we  saw. 

BOB  WEAVER  (UNIVERSITY  OF  SOUTHERN  CALIFORNIA):  What  was  the  total  extent  of  the 
temperature  fluctuations  in  your  graph?  I  could  not  quite  read  the  graph. 

MALCOLM  CALHOUN:  Minus  10  to  plus  35  -  yes,  minus  1 0  C  to  plus  35  C. 

OTOKAR  BUZEK:  These  peaks  are  fix>m  approximately  minus  18  to  plus  30  degrees  Celsius. 
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